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d In the context of quantum integrable systems (continuum or lattice), finding the non-Abelian infinite dimensional 

S [ algebra responsible of the existence of an (in)finite number of mutually commuting conserved quantities is a chal- 

ILJ ' lenge. Indeed, the success of conformal field theory essentially relies on the Virasoro algebra and its representation 

• i-H theory which allows to derive exact results such as the energy spectrum and exact correlation functions Qj . For 

KS lattice models such as the Ising 0, XY |2, superintegrable Z„— chiral Potts model 1^ and some generalizations 

j^ ■ |S], it is the Onsager non-Abelian Lie algebra |2] which plays an analogous role and allows to derive the exact 

spectrum for any conserved quantity of the model using solely its representation theory |Hj. In these latter models, 

the transfer matrix can be decomposed on a basis of mutually commuting conserved quantities /2fc+i, fc = 0, 1, ... 

which form an Abelian subalgebra of the Onsager algebra. In particular, each quantity l2k+i is expressed in 

terms of certain nonlinear combinations of the two fundamental operators ^Oi ^1 which satisfy the Dolan-Grady 

relations [5] 

[Ao,[Ao,[Ao,Ai]]]:^Po[Ao,Ai] and [Ai , [Ai , [Ai , Ao]]] ^ po[Ai , Aq] , (1) 

and generate the Onsager algebra with generators Ak,Gk, k integer [3 E] . Note that the parameter po as well 
as the realizations of Ao,Ai are model-dependent: for the Ising model po = 16 and for the superintegrable 
Z„— chiral Potts model po = v? . In both models, the generators are realized in terms of the SI2 loop algebra and 
the solution of the spectral problem for the transfer matrix or any conserved quantity has been solved using solely 
the representation theory of (^1 [HI . 

In the last few years, there has been some renewed interest in the Onsager algebra, Dolan-Grady integrable 
structure |H1 El 1101 1111 1121 [Qj and its deformation |14l 1151 [T?>j in the context of solvable lattice models as well as 
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1 Introduction 



in the representation theory of orthogonal polynomials P ll7[[TSlll9ll2Uj . Related with the subject of this paper, 
a g— deformed analogue of the Onsager algebra has been exhibited in various quantum integrable systems (Azbel- 
Hofstadter model, XXZ open spin chain with the most general boundary conditions ^t,...) which are usually 
studied using the quantum inverse scatterring method. As shown in |14II15| and similarly to the undeformed case 
12 El , in these models the integrability property is actually related with the existence of two operators Wo , Wi 
which satisfy the tridiagonal relations (also called the g— deformed Dolan-Grady relations) with scalars p, p* 

[Wo, [Wo, [Wo, Wi],],-i] = p[Wo, Wi] , [Wi, [Wi, [Wi, Wo],]g-i] - p*[Wi, Wo] . (2) 

More generally, they generate a g— deformed analogue of the Onsager algebra with generators {W_fe,Wfe+i, 
G/c+i, Gfe+i} introduced and studied in details in |lf)| . In addition, it has been discovered that the transfer matrix 
in these models - calculated in the inverse scattering framework |21j - can be simply written in the form 

t{u) = ^ T2k+i{u) l2k+i + J^oiu) I , (3) 

k=0 

where J-2k+i{u) are certain rational functions of the spectral parameter u and {l2fe+i}, fc = 0, 1, ..,N — 1, is a 
family of mutually commuting operators acting on the quantum space of the system. Remarkably, these operators 
form an ordered subset^ of the q— deformed analogue of the Dolan-Grady infinite hierarchy [21 El- Written in 
terms of the generators of the infinite dimensional g— deformed analogue of the Onsager algebra, this hierarchy 
takes the form: 

l2k+i -= K\N-k + K*Wfc+i + ^Gfe+i + ^Gfe+i , A: = 0, 1, ... (4) 

for arbitrary parameters k, k* , K±,k±. Although technically rather lengthy, by analogy with the undeformed case 
PI it is alternatively possible to write the operators {W_/j, W/j+i, G^+i, G^+i} for fc > in terms of nonlinear 
combinations of the fundamental ones {Wo, Wi}. For instance, one has ^21 (see also jTQ) 

Ii = kWo + K*\N, + ^[Wi, Wo]^ + ^[Wo, Wi]^ , (5) 

where the g— commutator [X, Y]q — q^^^XY — q^^/'^YX has been introduced. Note that the parameters iV, p, p* , 
K,K*, K±,k±, q and the functions J-2k+iiu) obviously depend on the physical characteristic of the model under 
consideration "^ . 

Having in mind the use of the Onsager algebra in the Ising '^, superintegrable Z„— chiral Potts model 2J or 
the Virasoro algebra in conformal field theory 1,, it seems to us important to analyse the models which enjoy the 
hidden dynamical symmetry |21l from the point of view of the representation theory associated with the tridiagonal 
algebra |2Il. Indeed, for several models - for instance the XXZ open spin chain with arbitrary integrable boundary 
conditions - the standard Bethe ansatz approach does not apply except at some special points in the boundary 
parameters space. As we are going to see, the analysis based on the representation theory of lO gives a clear 
understanding of the structure of the space of states in models with transfer matrices of the form (j^J. 

Among the interesting problems in integrable lattice models with hidden symmetry Q, one fundamental is 
to identify the eigenstates of the transfer matrix, and corresponding eigenvalues. According to the decomposition 
lO and the fact that the spectral parameter u is arbitrary, the structure of these eigenstates is encoded in the 
spectral problem for the g— deformed Dolan-Grady hierarchy {X2fc+i}, k = 0, 1,2, .... By virtue of |16| 

[I2k+i,l2i+i] ^ and X2fe+i1' = A2fc+i^ (6) 



^In the case of finite dimensional representations, the infinite hierarchy truncates due to the existence of linear relations among 
the generators. For the XXZ open spin chain, see 1161 . 

''The simplest example with a transfer matrix of the form lO corresponds to Af = 1: it is the Azbel-Hofstadter model which 
describes the problem of Bloch electrons in a magnetic field on a two-dimensional lattice. There, the elements Wq, Wi are realized in 
terms of the Weyl algebra ll4l . One has J^i (u) = v? , k, k* , k± , k± are related with length scales along the directions of the lattice and 
q ^ 1, related with the magnetic flux per plaquette, is a root of unity. An other model is the XXZ open spin chain with nondiagonal 
boundary conditions. In this case Wo,Wi are expressed in terms of U 1/2(^^2) operators acting on the Hilbert space of the whole 
spin chain with A'^ sites |14II15| . The parameters re, re*, re±, k± are related with integrable boundary conditions (see Section 3). 



for all fc, Z = 0, 1, ..., it follows that the structure of the eigenstates ^ is completely determined from the spectral 
problem for any operator 22fe+ii provided its eigenvalues are non-degenerate. For instance, one can focus on the 
simplest operator, namely Xi defined by |SJ|. If its spectrum is non-degenerate (see details in further sections), 
finding the eigenstates of Q is equivalent to solve Ii^ = Ai^. 

In this paper, we establish a direct relationship between the spectral problem for the hierarchy Q and a 
second-order g— difference equation in one variable - or its discrete version - satisfied by the states ^ in a suitable 
functional representation. The necessary and sufficient information to solve this spectral problem is shown to 
be encoded in the fundamental operators Wq, Wi, which properties uniquely determine the eigenfunctions '^(z). 
Some examples are considered in details. 

The paper is organized as follows. In Section 2, an infinite dimensional module of the tridiagonal algebra |(5J) 
is constructed. It is shown that the spectrum of the operator Wq can be algebraized i.e. corresponds to (possibly 
degenerate) polynomial eigenfunctions which roots satisfy a set of Bethe equations. In this eigenbasis of Wq, the 
action of Wi is described. Then, the spectral problem for © is considered: in the algebraic sector associated 
with special relations among the parameters n, k* , K±,k± ^ 0, the operator 2i admits polynomial eigenfunctions. 
Beyond the algebraic sector i.e. for arbitrary parameters k, k* ^ K±,k± ^ 0, the (non-polynomial) eigenfunctions 
^(z) of (0) are expanded on the (infinite dimensional) basis of polynomials eigenfunctions of Wq. Corresponding 
weight functions are determined by a coupled system of recurrence relations, which coefficients follow from the 
fundamental operator Wi. In Section 3, the same approach is applied to the case of finite dimensional modules of 
the tridiagonal algebra. In particular, we focus on the XXZ open spin chain with integrable boundary conditions. 
The eigenstates ^(-Zs) of the transfer matrix, represented by functions defined on a discrete support, are shown to 
satisfy a set of discrete q— difference equations. For general values of the boundary parameters, the eigenstates are 
constructed. For special linear relations among the boundary parameters, a subset of eigenstates (called Bethe 
eigenstates) are given by polynomial eigenfunctions defined on a discrete support, which roots are solutions of 
Bethe equations in agreement with known results. An alternative construction of these eigenstates is also proposed 
in Appendix A. Comments follow in the last Section. 

2 A second-order g— difference equation from the infinite hierarchy 

In this Section, we construct an infinite dimensional module of the tridiagonal algebra ^ on which the infinite 
g— deformed Dolan-Grady hierarchy Q acts. Although the infinite dimensional modules are more complicated 
to study than the finite dimensional ones, in our example the main properties of the finite dimensional case 
(eigenvalue sequences, block diagonal structure in the dual eigenbasis and duality of the fundamental operators) 
hold. These properties will be used to derive a g— difference equation for the eigenstates of Q. 

2.1 An infinite dimensional module of the tridiagonal algebra 

Let V denotes an irreducible infinite dimensional module on which Wq, Wi are both diagonalizable. Suppose q is 
not a root of unity. Let z denote a complex variable. Define r] to be the "shift" operator such that ri^^{z) = q^^z. 
Our aim is to endow the vector space of all functions in the variable x = z -\- z^^ with the module structure of 
a certain tridiagonal algebra |2Il. By analogy with the finite dimensional case (see ^JIEIEOI)) fo^' the linear 
transformations Wq, Wi it is natural to consider the following realization of the fundamental operators: 

Wo : cp{z)T^ + 0(z)7?-i + {d{z, z-^) - 0(z) - ^{z))I , 

Wi : z + z-i. (7) 

Here, <j){z)^ 4'{z), d{z, z~^) are functions of z which are not necessarely unique, as the spectrum of Wi may contain 
degeneracies. As we will see in Section 3, this kind of situation happens in the case of finite dimensional modules 
too. The existence of such degeneracies induces dim{V*) possible realizations of the operator Wq. 

*If <f>{z) 7^ 4'{z~^), we assume there exists at least one other reahzation of the operator Wq which spectrum is identical. 



with 



Our purpose is now to find the necessary and sufhcient conditions on (/)(z), 0(z), (i(z, z^^) such that the 
tridiagonal relations l(2Jl are satisfied. First, for any 4'{z), (t>{z), d{z, z^^) it is easy to check that Q automatically 
solves the second equation in ^ provided one identifies 

P* = -iq-q-')' ■ (8) 

Then, replacing {Tj) in the l.h.s. of the first equation of Q a straightforward calculation gives 

j{(j),^,d)(l)iz)irjix) ~ x)rj + r]-\-f{(j),^,d))'^iz){Tj-\x) - x)ri-^ (9) 

/?((/.,^,d) = {fj{q^{z))+r^(${z))-7j{d{z,z-'))){{q^-q)z + {q-'-q-')z-') 
+ (0(z) + ^(z) - d(z, z-l)) {{q - q-'){z-' - z)) 
+ (ry-i(^(z)) + v-'&i^)) - ri-\d{z, z-'))){{q-' - 9"')^ + {q - q^)z-') , (10) 

7(0,0,rf) = (77(<^(z))+r,(^(z))-r;(d(z,z-i)))' + (0(z) + 0(z)-rf(z,z-i))' 

-(9 + q-')[v{4>{z)) + viW)) - V{d{z, z-l))) (0(z) + 0(z) - d(z, z-l)) 

+ {2 + q + Q-i)0(z)7;(0(z)) + 0(z)r7-i(0(z)) + 7?(0(z))772(^(z)) 
, (l + g + g-i)/_ 



7y(x) — a; 



0(z)ry-i(0(z))(r;-i(a;) - x) + 77(0(z)),72(^(^))(^(^) _ ^2^^)^ ^ (^^) 



On the other hand. 



[Wo, Wi] = 0(z)(77(x) - x)77 + 0(z)(77-i(x) - a;)77-i . (12) 



Identifying ||5J) and (|12|) . we immediately deduce that the operators Wo,Wi with {Tj) satisfy the tridiagonal 
algebraic relations ^ iff the functions 0(z), (^(z), d{z, z"-^) solve the constraints: 

/3(0,0,d)=O and 7((?!),0,d) =p . (13) 

The equation (|10|) being a second-order g— difference equation with simple Laurent polynomials as coefficients, 
we restrict our attention to the solutions 0(z), (j>{z) which are rational functions of z and d{z, z~^) is supposed to 
be a Laurent polynomial in z. In this case, define the family 



2^_ ^ ^ and (/)(z) = </)(z ) , (14) 



-^ ' (i-^^)(i-9^^)n.=i (1-e.^) 

where we have introduced AN scalars {xk}, {Cfc} in ^^^ arbitrary field K. Replacing p4|) in 11U|) . (|ll|l . eqs. (|13|l 
impose some constraints on the parameters: 

For A^ = 1, it is straightforward to check that the constraints p3|) are satisfied with the identification ^ and 

p^ -{q- q^^fxiX2X3X4q^^ , d{z, z^^) = i + xiX2X3X4q^^ , (15) 

for any choice of the parameters {xiiX2,X3jX4}- Notice that the realization Q for the operator Wq in this case 
coincides with the second-order Askey- Wilson g— difference operator |22| . 

For N = 2, the first constraint in (|13() is satisfied for the choice 

d(z, z"^) = 1 + XiX2X3XiX5X6Ci^^2^q~^ , (16) 



whereas the second constraint yields to the foUowing relations among the parameters: 

3 2+fc 

fc=-2 ii<...<i3+fc j=0 

E( E (-i)'=xnx....x.3..(a6r'(Eci6^^'^^)) = 0, (17) 

k= — 3 ii<...<i3+k j=0 

with {ifc} G {1, ...,6}. So, it remains six independent parameters, say Xii X2, X3: X4: X5: Xe- Simplifying (|13|l . one 
finds 

p = -{q- q~^fxiX2X3X4X5X6S.T^^2^q~^ ■ (18) 

Although we do not proceed further'^, for higher values of N it is possible to classify all possible parameter 
sequences such that {Tj) with p4|) and d{z,z^^) = Const, satisfy jSJl. More generally, it is possible to find 
examples for ip, (p which are neither such that (j>{z) = (/)(z^^), nor rational functions of z. 

2.2 A polynomial eigenbasis 

For the rational functions 0(z),0(z) satisfying (|13() . a basis of possibly degenerate eigenfunctions ipn[m]{^)i '^ — 
0,1,2,... and m — 1 , . . . ^ dimiVn) of Wq with eigenvalues {A„} can be constructed. The realization (|7|l for Wg 
yields to the g— Sturnr-Liouville problem 

(t>{z)i),,[,n] iqz) + '4>{z)i-'n[m.] {q'^ z) + fi{z, 2:"^)V'„[m] {z) = A„ Vn[m] (z) (19) 

with 

/i(z, z~ ) = d{z, z^ ) — (f>{z) — 4>{z) . 

In general, only a part of the spectrum of (7— difference operators can be algebraized, i.e. corresponds to 
polynomial eigenfunctions in z. However, according to the structure of the spectrum of Wo;Wi for the finite 
dimensional case ^|E1) we expect the operator Wq in Q can be entirely algebraized provided (fU^ are satisfied. 
For this reason, we are looking for polynomial eigenfunctions symmetric in the variable z i.e. invariant under the 
replacement z -^ z~^: 

n 

i^nlm]{z) = l[{z-z,){z-^-z,) for z, gIzJ^^}, (20) 

with m — 1, ...,dim{Vn), n = 0,1,2,... and Zj denote the roots of the polynomials. Replacing this expression 
in ifliHl and setting z ^ Zi, a,n immediate consequence of the factorized structure ipn|l is the system of Bethe 
equations 

^"('^"^) - ^("^) with . = l,2,...,n. (21) 



^n{q ^Zi) 0(Zj) 

As an example, let us consider the family of rational functions H14|l with suitable parameter sequences such 
that H13|) are satisfied. Introducing the new parametrization 



z = e^\ XI = e^"' , 6 = e'^' , q = e^ , (22) 
the Bethe equations take the form 

-r-r sinh(A^ + q) -p-r sinh(Ai - r/i) _ -A- sinh(A.t + Xj + ip/2) sinh(A., - Aj + ip/2) 

fi sinh(A, --q) fjl smhiXi+r]i)~ _\^ smh{X, + Xj^ip/2)sinh{Xi^Xj-ip/2) 



^The case N = 3 has been also considered, although not reported here. 



with i — 1,2, ...,n. For A^ = 1 and no special relations among the parameters {xk}, one has dim{Vn) = 1- The 
polynomial eigenfunctions H2U|) coincide exactly with the Askey- Wilson polynomials which zeros are known to 
satisfy H23|) at A^ = 1 . For N — 2 the parameters ci , C2 are determined by the relations H17|) , and the solutions 
1)20(1 can be seen as some generalization of the Askey- Wilson polynomials which, to the best of our knowledge, are 
new. In both exactly solvable cases, the spectrum can be easily derived. Replacing ()20|l in H19|l. expanding both 
sides of the equation and identifying the leading terms one finds 

K = xiX2X3X4g"^9" + g"" for N = i , 

K - XiX2X3X4X5X6Cr'C2"''?''9" + '?"" for N = 2. (24) 

For higher values of TV, a spectrum of the form 

y^0i^nlm]{z) = Xntpn[m]{z) with A„ = Cg" + C'^"" , (25) 

C = 9{{Xk},{ik},q), C" = 1 is clearly expected: the eigenvalue sequences for iV = 1,2 and N = 3 although 
not reported here agree with H25|l . More generally, such form of the spectrum coincides exactly with the one 
associated with finite dimensional modules |17lll8j (see also |2()j'l. 

Let us now consider the action of Wi on the eigenbasis of Wq. According to iQ and the fact that f/'nfml (-2^) 
are symmetric polynomials of degree n in x — z + z^^ , the operator Wi has a block tridiagonal structure in the 
eigenbasis of Wq. Its action on ipnlm]{z) can be written formally® 

dim{Vn^i) dirn{Vn) dim{Vn-i) 

WiV'„[m](^;) = ^ b.n[l„i]i'n+l[l]{z) + ^ Cl„[hn]i^n[l]{z) + ^ Cn[im]ll'n-l[l]{z) ■ (26) 

1=1 1=1 1 = 1 

For the family of rational functions (|14l) and d{z,z^^) = Const., the coefficients Cin[hn]-> OnUm]} Cn[im] can 
be determined explicitly in terms of the parameters Xfc,^/c- As the simplest example, consider the fundamental 
operators Wo,Wi with and !(T^ for A^ = 1. To avoid degenerate situations, we assume the parameters {xfc} 
are generic. As mentionned above, in this case dim(Vn) — 1 i.e. m,l = 1. The coefficients in ((26|l are well-known 
|22| . they take the form 



^n 



111 



(1 - xiX2g")(i ~ xiX3g")(i - xiX4g")(i - xiX2X3X4g" ^) 

xi(i - xiX2X3X4g^""^)(i - xiX2X3X4g^") 
xi(i - g")(i - X2X39"-^)(i - X2X4g"^^)(i - X3X4g"-i) 



(1 - XlX2X3X4g^"' ^)(l -XlX2X3X4<?^" ^) 

O-nlu] = Xi + xr ^ L>„[ii] - C„[ii] . (27) 

To resume, polynomial eigenfunctions of the form ((2()|l with Bethe equations H21|) provide an example of 
irreducible^ infinite dimensional module of the tridiagonal algebra. In this basis, Wg, Wi act as (|25|l and H26|) . 
respectively, where explicit expressions of the coefficients depend on the solutions (j){z), (j){z),d{z, z"^) of (fO|l . For 
the family of rational functions H14|l , the Bethe equations take the form l(23|l with the identifications ((22(1 . 

2.3 Structure of the eigenfunctions 

For quantum integrable models with dynamical symmetry ((SJ, the spectral problem for the family of mutually 
commuting operators {l2k+i} k = 0, 1, ..., reduces to the one for (0) provided the spectrum of this operator is non- 
degenerate. Let ^(2) denotes an eigenfunction of the hierarchy Q. Using the realization jT)) of the fundamental 
operators, the spectral problem for Ii in the functional representation reads ^ 

A{z)^{qz) + 'A{z)^{q-^z) + B{z, z"^)*(z) = Ai*(z) (28) 



^This property can be also derived starting from J2J. See for instance 1181 . 

^For special relations among the parameters entering in the solutions (t>{z), (j>{z), d{z, z~^), the module may become indecomposable. 
This possibility is not considered here. 

*The reader familiar with the algebraic Bethe ansatz will immediatly recognize the similarity of this equation with the Baxter 



identity for the eigenvalues of the transfer matrix, also called T-Q relations [22 



with 

A{Z) ^ («+(g-,-l)(^g-V2,-l + ^^l/2,))^(,)^ 

V fc+ A;_ / 

Here, the parameters k, k*, k±, A:± are arbitrary. Contrary to the case of Wq, the operator Zi admits only partial 
algebraization of it spectrum. Depending on the parameters, we now describe the structure of its eigenfunctions 
\E'(z) in the non-algebraic sector and the algebraic sector, respectively. 

• Non-algebraic sector: For generic values of the parameters k, k* ,K±,k±, there are no polynomial solutions. 
Indeed, the leading contribution in the numerator of l|28() - independent of Ai - can not be cancelled out. A natural 
approach is then to consider ^(z) as an infinite power series expansion in the variable x = z + z^^. However, 
for the family of rational functions p4|l the problem becomes quickly rather complicated when the value of N 
increases. Instead, a convenient procedure consists in expanding the solutions of (|28|l on the basis of eigenfunctions 
4'n[m]{z) of Wq givcn by H2UII with ]2'A\i . For generic values of the parameter k, k*, k±, k± we write 

oo dini{Vn) 

'^(^) = Z1 Yl fnlrn]{^l)i'n[m]{z) ■ (29) 

n— m— 1 

According to (0 and the action H25(l of Wq and l|26|l of Wi on the polynomial basis, the coefficients /„[m](Ai) are 
determined by the following coupled system of recurrence relations: 

dini{Vn-i) dini(Vn) dim{Vn + i) 

[nXn - Ai) fn[l] + 2^ Sn-l[lm]fn-l[m] + 2^ Ai[im]/n[m]+ 2^ C„+ 1 [;„]/„+ 1 [„]= , (30) 

771 — 1 777 — 1 777—1 

for / — 1, ..., dim{Vn), where the coefficients Aniim], l3n[i,n],Cnihn] follow from the action of © on (|2()|l . Using (|25|l 
and (|2SJ), they read 

^771.777] = {^* + i^q'^' i^'^"'^')^" + ir'''^' |7'^"'^')^" + l)) ^"['™1 ' 

C77[,7,7] = {^* + {^q'" ^<i-"')^n + {'^q'l' ^q-'/^)Xr..,)) C„[,,„j , 

A.n[i,n] - ('^* + (g'/'-g-'/')(^ + ^)A77)a77[i777] . (31) 

• Algebraic sector: This sector corresponds to polynomial eigenfunctions of Xi which exist for special 
relations among the parameters, identified as follows. Let us assume that ^71(2:) is a polynomial symmetric in the 
variable z of the form 

77 

*77(^)=n(^-^^-)(^"'-^.), (32) 

with n = 0, 1, .... For the family of rational functions p4|l with suitable non- vanishing parameters, replacing (|32|l 
in H28|) one finds that the leading contribution in the l.h.s. of the equation (obtained by taking the limit z — > 00) 
is independent of Ai. Then, the leading terms will cancel out (and similarly if one considers the limit z —f 0) only 
if the parameters satisfy the relations 

2JV+2 2N-2 

^* + iq-q-')^q-'^'q-'' + {q-q~')^q'^\q-' [I ^^ 11 ^k)<l" = ^^ (33) 



i.e. Bn[im] = 0. This implies that any eigenfunction H32|) can be written as a finite sum of elementary eigenfunctions 
H2U|) . Choosing for simplicity the representation ()32|l . the roots Zi are now determined by the generalized system 
of Bethe equations 

— tor 1 — 1,2, ...,n . (34) 



Given n fixed by (|33|l . the spectrum can be calculated by plugging (|32|l into (|28|l and identifying the leading terms 
of both sides of the g— difference equation. For the family of rational functions H14(l . the result is 



i=l 



,{N) AN) 



where T^ iQ± are some functions of the parameters {xk\^ {^fe}- For small values of N , 

ii<i2<i3 i=l 

6 

H'^ = XlX2X3X4X5X6er'e2"'9-' , ^i'^ = 1 , ei'' = 5] X^ " Cl ^ 6 , 

4 

^+' = ( Yl XtiXt2Xt3Xt4Xt5)^1^^2^ -XlX2X3XiX5X6^1^^2^l'^\^l^ +^2^) ■ 

il <i2 <i3 <i4,<i5 

For higher values of N, a similar form of the spectrum is expected according to (|14|l . It is however important to 
remind that the parameters {^k } are restricted by some relations generalizing p7|) . For non- vanishing parameters 
K±, the spectrum H35|) is non-degenerate. So, (|32(l is also a polynomial eigenfunction of the infinite hierarchy Q. 

Consequently, the operator Xi admits partial algebraization of its spectrum. In the algebraic sector associated 
with (|33|l . there is a single eigenstate associated with the polynomial eigenfunction l|32|) with 134|1 . For the family 
of rational functions p4|l . the spectrum can be written in terms of the roots of the Bcthc equations, and takes 
the form (|35l) . The rest of the spectrum is associated with non-polynomial eigenfunctions which can be written 
in the (infinite dimensional) eigenbasis H2U|) with \21\i of Wq . The structure of the eigenstates (|29f) is encoded in 
the coefficients /„[,„] which satisfy H3U|) . 

3 A discrete g— difference equation for the XXZ open spin chain 

For quantum integrable lattice models with dynamical (7— Onsager symmetry, a similar analysis can be done in 
order to derive a g— difference equation for the eigenstates of the transfer matrix Q. The main difference is 
however the Hilbert space of these models, which is now finite dimensional^. To give an illustration, let us 
consider the XXZ open spin— 2 chain with general integrable boundary conditions. Its Hamiltonian reads 

U S^ ( „k„k+l I k„k+l I \„k„k+l\ 

tlxXZ = 2^ (^O-jCTi +0-2 0-2 +^Cr3CT3 J 

fc=l 

+ K^l'T^ (^^^a" + 2(.^/^ + <l-'''){n,a- + ._ajV)) , (36) 



^Linear relations among the higher operators of the g— deformed analogue of the Onsager algebra - generated from Wq, Wi - exist 
that are responsible of the truncation of the hierarchy Jl) (see |16| for details). 



where cri.2,3 and (7± = (cti ±i(T2)/2 are usual Pauli matrices. Here, A — {q^^^ +q ^^^)/2 denotes the anisotropy 
parameter with q = exp (p. Integrable boundary conditions correspond to the choice j24| 

ef = (eW)t ^ (coo + iCoi)/2 , k+ - (fc_)t = -(gi/2 „ q-i/2y0/2 (left) , 

K = {k^ = (-5oo + *5oi)/2 , «+ = (a.-)^ - -e'V(2(<zi/2 + ^-1/2)) (^gi^^) . (37) 

The r.h.s parametrization has been introduced in which case one immediatly identifies the Hamiltonian as defined 
in |25| . In total, one has sia; boundary parameters Cqo,Coi,Coq,Cqi,6,0. Except if explicitely specified, from now 
on we assume these parameters are generic. 

As discovered in |15[ [TU| , the known integrability of the XXZ open spin chain (|36|l with (|37|l |21l ^^ is related 
with the existence of a hidden dynamical symmetry, a g— deformed analogue of the Onsager algebra. This non- 
Abelian algebra is generated by a family of nonlocal operators {V^^k ,^k^i^^k+i^^k+i}^ which fundamental 
ones satisfy the tridiagonal relations ^. They read JSj 

Wi""^ = (fc+a+ + fc_a„)®J(^-l)+g-'^■^/2 0>vf-'^ (38) 

with "initial" conditions 

W^''> = fc+a+ + fc_a_ + ef g'^-^/2 , 

W['^ - k+a+ + fc_<7_ + eL"^g-"-^/' . 

Note that the left boundary parameters in (|35|l appear in the realization (|38|) of the algebra Q, whereas the 
right boundary parameters enter in the definition of the linear combinations Q. For the normalization above, 
the structure constant of the tridiagonal algebra |5J is fixed to 

P = P = 1 • (39) 

Irreducible finite dimensional modules V of the tridiagonal algebra |(2Jl associated with Wq , W} have been 
studied in details in (SO]. Here, we just recall some results. Suppose the boundary parameters entering in (p?H|l 
are such that both operators are diagonalizable on V, and suppose q = e'^ - with ip purely imaginary - is not a 
root of unity. Then, according to [JHI) Theorem 3.10] the operators Wq , W^ act on F as a tridiagonal pair: 
The pair of linear transformations Wq ' : V ^ V and WJ; . V ^> V satisfy the following 

(i) There exists an ordering Vq,Vi, . . . , Vn of the eigenspaces of Wg such that 

>VrV„ C K+i + K + K-i (0<n<iV), and F_i = 0, Vat+i = , 

(m) There exists an ordering Vi^ ,V* , . . . , V^ of the eigenspaces of W^ such that 

W^^V; C f;+i + V: + VU (O < s < iV), and T/_*i = 0, 1/^+1 = . 

Also, for < n < A'^, the dimensions of the eigenspaces Vn and V* are equal and the shape vector '\T of the pair 
is symmetric and unimodal. For the family of operators H38|l , one has ^° |3Dj 

dim{Vn)=y\ Q<n<N . (40) 

Depending on the basis chosen the matrix representing Wg (resp. W} ) is either diagonal (with possible 
degeneracies) or (block) tridiagonal in the dual basis. Below, we will use these properties and refer the reader to 
the work j2()| for the explicit expressions of the matrix representation in each basis. 



^We denote 



JV! 
n!(iV-n)! 



3.1 Functional representation 

For our purpose, it will be useful to consider the functional representation of dimension dim{V) — 2^ defined 
on a discrete support associated with the overlap coefficients between the two dual eigenbasis of Wq , Wf . 
As follows from (i) and {ii), these coefficients satisfy a set of coupled recurrence relations and a set of discrete 
g— difference equations, respectively. Introducing the parametrization 



e^ — cosh a and e_ — cosh a* , (41) 



they are denoted 

V'S (^-) ^itli ^s - exp(a* + (TV - 2s)ip/2) , (42) 

where n = 0,1, ...,N and ?7i = 1, ..., ( ^ ) characterizes the degeneracy of the eigenvalue with index n. Note that 
these coefficients are symmetric in the variable Zs, as we will see below. They provide a functional eigenbasis for 
W^^\ Indeed, one has EHI 

WrVJ™](^«) = ^nV'i7i(^^) ^^*^ A„ = cosh(a + (TV - 2n)^/2) . (43) 

On the other hand, the action of the operator W]; in the basis H42|l can be explicitely calculated. It takes a 
block tridiagonal structure which is formally written as 

1=1 1=1 1=1 



where the explicit expressions of the matrix entries can be found in |2()j . Analogous results also hold if, instead, 
we had considered the dual overlap coefficients which also form a complete eigenbasis of WJ , a phenomena due 
to the duality property of (01. Obviously, in this dual basis the operator Wq has a block tridiagonal structure. 

The analysis of previous Section can be applied to the finite dimensional case similarly: the operator Wq 
can be associated with a discrete g— difference operator acting on the (discrete) functional space spanned by 142|l . 
Defining iT^^{f{zs)) = f{zs±i) = f{q^^Zs), the discrete analogue of the realization O is now given by (^ 



Wr : <'\s)Vs^^[''\s),:^+,f\s) , k e {!,..., (1 



Wj^) : (z, + z-')/2 , (45) 



where 



1=1 ^fc I'^J 

The explicit expressions of the coefficients can be found in 20 . It is important to notice that the realization of 
Wq is not unique (except for s = and s = N), as the spectrum of W} contains ( ) degeneracies. 



3.2 Eigenstates of the XXZ open spin chain 

Due to the non-Abelian q— Onsager dynamical symmetry of XXZ open spin chain (|36|l . the transfer matrix can be 
written in the form (j3l with Q J15II16| . By analogy with the case of infinite dimensional modules, all necessary 
and sufficient information to determine the eigenstates of (jSJ is encoded in the structure of the tridiagonal pair 
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^0 1^1 given by Ip!^ . We choose to take l|l^ as the basis on which any eigenstate, denoted ^(zs), is 
expanded. Considering the spectral problem for Zi and using the realization (|45|l . we find that the eigenstates 
of the XXZ open spin chain with general boundary conditions (|36|l are the family of solutions to the following 
system of discrete g— difference equations 

Af'(z.)vI/(9z,)+4'^^(z,)vI/(g-iz,) + sf)(z„z7i)vI/(z,)=AiVl/(z,) , fc-l,...,(^) (46) 



with 



Bf\zs,z-') = -A^f'(^)+(«* + (.^/^-.-/^)(^ + ^)Mr(^))^^^^, 

gni and s = 0,l,...,iV. 

Difference equations are known to be closely related with discrete equations. For the example considered here, 
the q— difference equation ^^ (|28|l can be seen as an extension of (|4mi on the whole complex plane, where the 
parameters {xfc}, {S,k} in (|14|l are chosen such that 

0(z = z;i) ^ ^[''\s) , 0(z = z-')^^f\s) , M(^,^-')L=.r^ - ^^(5) , -^{z = z:') ^ ^{zs) . (47) 

Indeed, finite dimensional module of the tridiagonal algebra (0) can be obtained from the infinite dimensional 
modules by imposing z to be on a discrete support z — Zs, s = 0, 1, ..., iV. It follows that the results of previous 
Section can be used to analyse the structure of the eigenstates of the model (|36|) . In particular, in view of H33() . one 
expects that these eigenstates admit a description in terms of solutions of Bethe equations for certain relations 
between the boundary parameters. Below, we call the corresponding special cases the algebraic sector of the 
operator (01. 

• Non-algebraic sector: For generic values of the boundary parameters k, k* , k± , A;± , there is no non-trivial 
subspace of V which is left invariant by the action of ©. So, any eigenstate is a linear combination of all the 
eigenfunctions H42|) . n = 0, 1, ..., A^. In general, it takes the form 

*(^«) = E E /nM(Ai)<i](^.) ■ (48) 

71—0 Tn — 1 

The eigenfunctions (|42|l being clearly identified |2n|, it remains to determine the non- vanishing coefficients fn[m]- 
Using the block tridiagonal structure (|44|l of WJ in the functional basis H42|) . we find that the coefficients 
/n[m](Ai) satisfy the set of coupled recurrence relations (|3(J|I with the definitions H43|) . (|40|) and the explicit 
expression for the coefficients taken from jJOj. Note that each eigenfunction (|20|l at z = Zs is individually 
associated with a system of Bethe equations, whereas there is no single system of Bethe equations for ^(zs). 

For small values of N , it is not difficult to check the recurrence relations %M)\ numerically ^^ . One finds that the 
2^ values of Ai agree exactly with the direct diagonalization of Xi for generic values of the boundary parameters 
H37(l . As expected, the corresponding eigenvectors diagonalize the higher operators Q. 

• Algebraic sectors: The advantage of the formulation of the transfer matrix in terms of Q) is that the 
representation theory of the tridiagonal algebra Q guarantees that the discrete g— difference equations (|46ll admit 
a subset of solutions associated with Bethe equations. This subset is identified with the special points in H28|l with 

-'^-'^The class of g— difference equations I28i may be extended by a "gauge" transformation ^(z) -^ f{z)'if{z), (j>{z) -^ (f){z)f{z)/f{qz), 
<f>{z) -^ <f>{z)f{z)/f{q~^z) and fi{z) -^ IJ-{z), where f{z) is a rational function of z. 
^^K. Koizumi, private communication. 
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103); Pn|) - or equivalently l|l^ using (|17|) - where either some of the coefficients A)^ (zs), A,, (zs) or Bn[m]j ^^[m] 
in (|31|l are vanishing, or the special points for which the eigenfunctions (|42|l are no longer linearly independent. 

For instance, suppose Bn[im\ = 0. Plugging the parametrization (|37() . H41() . (|43|l into (|31|l . one obtains a linear 
relation between the left and right boundary parameters which reads 

a + a^i{e~e)-{N -2n+l)ip/2 mod{2iTi) . (49) 

According to the action of the operator Zi on the eigenfunctions of Wq : 

Vp+i; \p) Vp-1/ 

;=i i=i 1=1 

one has 

Then, 

IiWn C W^„ where M/„ = F^ 

p=0 

is an invariant eigenspace of Q. The eigenstates ^^^((zs) associated with H32|l at z = z^ and (|34|l form an eigenbasis 
of Wn which dimension follows from the degeneracies of each eigenspace of Wg , given by 1)40(1 . One has 

d*m(W„) = Vf^V (51) 

This invariant subspace of Q) is not unique: another invariant subspace 

N N 



Wn^^Vp with dim{Wn) = X ( ^ 

p— n p— n ^ 



exists provided the boundary parameters are suitably tuned such that Cn[im\ = 0. Using (|^ . |(1T|) . ll^ in l(^ . 
this condition corresponds to 

a + a = i(6'-^) - (Ar-2n- l)v7/2 mod{2iTi) . (52) 

In this case, using the substitution n — * iV — n in (|32|l at z = Zs, the eigenstates with ()34|l form a basis of 

Wn- It is important to notice that the conditions A^, (z^) = or Aj^ [zg) = may have been considered 
alternatively, in which case the linear relations between the boundary parameters are obtained from (|49() , H52() by 
complex conjugation. This is not surprising, in view of the duality between the eigenbasis of (|38|l . To complete 
the analysis of the eigenstates in both algebraic sectors 1(491) . 1(521) . it is worth mentionning that the solutions 
^„(zs) can be regarded as linear combinations of the eigenfunctions ((4211 . truncated in comparison with ((48() . The 
coefficients in the expansion satisfy a set of recursion relations, reported in Appendix A. 

Other invariant subspaces exist if the eigenfunctions ((42(1 entering in ((5()() are not linearly independent. These 
are the special points where the spectrum ((43(1 - or the one in the dual basis - admits additional degeneracies, 
namely 

a + {N - 2n)Lp/2 = ^ttZ , (53) 

in perfect agreement with |26[ . The same result can be also derived from the explicit form of the eigenfunctions 

^n[m]{Zs) M- 

To conclude this Section, the approach based on the tridiagonal algebra (O provides a new derivation of the 
linear relations ((49(1 . I(52() . ((53() . In our framework, the (Bethe) eigenstates "^nizs) correspond to the algebraic 
sectors of the q— deformed Dolan-Grady hierarchy Q), where the Bethe anstaz equations have been formulated 
and the spectrum of the XXZ open spin chain has been derived [23 1251 |2H| . 
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4 Comments 

A new approach based on the algebraic properties and irreducible modules of the infinite dimensional g— deformed 
analogue of the Onsager algebra has been proposed in the context of quantum integrable models. The main 
result of this paper is an exact second-order g— difference equation (resp. its degenerate discrete version) for the 
eigenfunctions associated with the infinite (resp. truncated) g— deformed Dolan-Grady hierarchy Q. For generic 
values of the parameters k, k* ,k±, k± and variable z (resp. discrete Zg), the structure of the eigenfunctions has been 
described in details. They admit an infinite (resp. finite) serie expansion in terms of the polynomial eigenfunctions 
of Wq which roots satisfy a system of Bethe equations of the form (|23|l . whereas the weight functions satisfy a 
coupled system of recurrence relations which explicit form is determined from Wi . In the algebraic sectors which 
correspond to special relations between the parameters, the expansion truncates. In this case, the eigenstates 
admit a formulation in terms of solutions of the system of Bethe equations H34|l . The approach presented here has 
been applied to the XXZ open spin chain, although higher spin representations can be obviously studied similarly. 

Equations of the type (|28() . (|46|l already have important physical applications. For instance, the method of 
algebraization has been applied to the Azbel-Hofstader problem ,29] (see also PU])- In this model, the spectral 
problem for the Hamiltonian (Schrodinger equation) H = 2i with anisotropy parameters k, k* is the famous 
Harper's equation. For certain values of the physical parameters, it can be algebraized using its connection with 
Uq{sl2) at root of unity. Non-polynomial solutions have been considered in j31| . 

Such a connection between g— Sturm-Liouville problems and quantum integrable systems is not new. However, 
its understanding and formulation using the g— deformed Onsager algebra/tridiagonal algebra (j^J studied here is 
rather promising. A classical counterpart of our description clearly needs further investigations. In this direction, 
let us mention that a classical analogue of (0) has been studied in 32 . The simplest example of realization of 
classical Leonard pairs Wq, Wi has been proposed: 

Wo = y^ + U{x) and Wi == 0(a;) with {x, y} = 1 , (54) 

where the admitted potentials U{x) are the (i) hyperbolic, (ii) modified hyperbolic {Hi) trigonometric Poschl- 
Teller potential, or the {iv) Morse, (v) singular and (vi) shifted oscillator potential. Indeed, in the continuum 
limit g ^ 1 of l(3) with l|14|l at iV = 1 the g— difference equation H28|l turns into the Schrodinger equation for the 
generalized Poschl- Teller potential. Furthermore, using a suitable change of variable and "gauge" transformation 
one arrives at the Heun equation. In view of the recent works j33| . we expect our approach will provide an 
algebraic understanding of the ordinary-differential equations/integrable models (ODE/IM) (for a review, see [51] 
and references therein). In the quantum case, it should be also mentionned that our results exhibit some links 
with a recent work relating a g— Sturm-Liouville problem with the Bethe ansatz equations of the XXZ open spin 
chain with Dirichlet boundary conditions |12| . In particular, for vanishing parameters {^k} the rational functions 
(|14l) considered here can be easily related with the ones in |3Hj. However, for this choice the g— Onsager dynamical 
symmetry is broken - as expected - and {T)) do no longer satisfy ^. 

To conclude, the program initiated in 14; opens the possibility of studying massive quantum integrable models 
from an algebraic point of view (model-independent) which idea takes its roots in the original work of Onsager 
12]. The exact spectrum of the complete hierarchy Q) is the subject of a separate work, among other problems 
mentionned above. 

Note added: The explicit expressions for the entries of the matrices ()38|l in their dual eigenbasis are not 
reported here. We refer the reader to [201 for these data. 

Acknowledgements: I thank P. Forgacs, H. Giacomini, N. Kitanine, K. Koizumi and R. Weston for discussions, 
and P. Terwilliger for suggestions and interest in this work. I wish to thank the organizers of the 3rd Annual 
EUCLID Meeting 2005 where preliminary results were presented. Part of this work is supported by the ANR 
research project "Boundary integrable models: algebraic structures and correlation functions^'' , contract number 
JC05-52749 and TMR Network EUCLID "Integrable models and applications: from strings to condensed matter'^ , 
contract number HPRN-CT-2002-00325. 
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Appendix A: An alternative construction of the Bethe eigenstates 



For special relations between the left and right boundary parameters, the corresponding Bethe eigenstates can 

eigenfu 

n \ p J 



be written directly in terms of the eigenfunctions of Wg '' . Consider the linear combination 



*n(^«) = EE/pH(Ai)V4S(^^) ^0=^ n<N-l, (55) 



where /'r^i are non- vanishing coefficients. As an immediate consequence of H5U|) . given n this combination is an 
eigenfunction of Xi iff the condition Bniim] = i.e. (|49|l and the coefficients for p = 0, f , ..., n — f with degeneracies 
; = 1,...,(^) satisfy 

Vp-1/ \p) Vp+1/ 

m—1 772—1 m — 1 

and 

(kA„ -Ai) /,;[;]+ ^ ^„[;„]/;[„] + Y^ C^+lllr,i]fn+llrn] ^ ^ for ^ " -^' - '' ( „ J ' ''^''^ 

771 — 1 777 — 1 

Antoher choice is Cn[im] = Oj which corresponds to the linear relation between the left and right boundary 
parameters H52f) . An eigenfunction of Xi takes the form 

N (p) 

^»(^^)-EE/pH(Ai)^iM(^«) fo"^ "^1' (58) 

P— 77 777—1 

where the coefficients in the expansion for p ^ n+ 1, ..., N with degeneracies Z = 1, ..., ( ) are determined by the 
set of recurrence relations 

Vp-ij \p) Vp+ij 

[kXp - Ai)/pj;] + 2_^ Bp^i[hn]fp-i[rn] + E "^pCH /p[7n] + E '^P+lfM /p+lM = , (59) 

777 — 1 777 — 1 777 — 1 

and 

(kA„ - Ai)/;^,] + Y^ I3n-llhn]fn-l[m] + E A7 [/777] //j [777] = for I = 1, ..., ( j . (60) 

m^l 777—1 



Appendix B: Identical polynomial eigenfunctions for the descendents 

The purpose of this Appendix is to give further support on the fact that higher operators of the form W_fe, 
for instance W_i, are diagonalized by the polynomial eigenfunctions (|2()(l . in agreement with [Wq, W_fc] = for 
all A: = 0,1,.... One has [El 

W_i = -i[Wo, [Wo,Wi]J,-i+Wi and W2 - -^[Wi, [Wi, Wo];,-i + Wq . (61) 

Using previous results, suppose {Wo,Wi} take the form {Tj) and assume 0,0 and d solve (fT^Hl . Replacing in 
1)61(1 . after simplifications one finds the following action on the infinite dimensional module 

W_i : 0_i(z)77 + 0-i(z)?7-i + //-i(z,z-i)/ , 

W2: i^2iz,z-^) (62) 
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with 

-1\2 



0_i(z) = -i(((q - g2)2 + (g-i _ q-^)z-^)^,{qz, q-^z'^) + {{q-^ - l)z + (1 - q)z-^)li{z, z-^))<P{z) , 
^_i(z) = --(((q-^ - q-^)z + {q- q^)z-') fiiq-' z, qz'^) + ((1 - q)z + (1 - g-l)z-l)Ai(z, z-^))'4>{z) , 
/i_i(z, z-l) = --(((92 - l)z + (g-2 _ l)z-l)</,(z)^(qz) + ((g-2 - l)z + (g2 _ l)z-i)^(z)^(g-iz) 
+ (z + z-i)((g + g-i - 2){^l(z, z-^)f^ + z + z'K 



The tridiagonal structure H62(l is not surprising: similarly to the finite dimensional case"'^'^, the action of higher 
order operators can be written solely in terms of the g— difference operators 77 and the identity /. The g— Sturm- 
Liouville problem associated with W_i takes the form 

0_i(z)^(qz) + 0_i(z)V(g"'2) + /i-i(^, z^^Mz) = AV(z) (63) 

with the definitions above. Polynomial solutions can be constructed similarly to the case of Wq. It is however 
important to notice that, for any z, 

^-i(^) = M. , (64) 

(j)-i{z) (f){z) 

an immediate consequence of HlOfl and l|13(l . Then, H64|l implies that the system of Bethe equations associated 
with Wo is identical to the one for W_i. Consequently, up to some overall factors polynomial eigenfunctions of 
both operators coincide. 
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